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IJNTIiOnUCTION 


Tlie correct selection of dynamic characteristics of n nuclear power plant will allow to 
provide its maximum reliability in the transient and accidental conditions. 

One of the most essential branches in the dyiininics of reactor systems of a water-moderated 
water-cooled type is the analysis of transient processes at a chanf^ of the coolant circulation 
in the primary circuit. In this event the followiiu. operaiint^ conditions of the |)Ower plant are 
of a great practical interest: 

(a) the accidental drop of a coolant flow rate; 

(b) the increase of a coolant flow rate at the pump starting. 

The drop of a coolant flow rate may result in an unnecessary abrupt overheating of the fuel 
element cannings even after the chain reaction is ceased. Ivspecially dangerous is melting of 
the fuel element canning due to the possibility of fission-product release into the circuit. 
Besides, at the canning tempemturos somewhat higher than the melting point, the melting 
metal dangerously interacts with the water or steam generating a great amount of 
heat [!]• 

When starting the pump some sharp uncontrolled surges of neutron (and heat) power of the 
reactor may occur at a usually negative temperature coefficient of reactivity. 

In this case the power increase is connected with the release of positive reactivity at a 
decrease of the bulk temperature of water in the reactor (the decrease of water temperature in 
the reactor is caused both by the increase of a coolant flow rate and by the possible delivery 
of some colder water from the starting pump pipeline to the reactor). 

The designer's problem is to define the correlation between the characteristic parameters of 
a process leading to the predetermined transient conditions. The predetermined transient 
process may be characterized by a tolerable enthalpy of the coolant in the reactor and by the 
lack of heat exchange crisis (or by the maximum tolerable temperature of the fuel element 
canning). 
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MathemnUcally tha problem \h eonfinod U> « Hyatem of .lifferonlial e(|u.ilion« of common 
variable ilerivalivcH winch an- (|.artially) obtaiuud by findinK »!•« mean values of equations in the 
purtitti dorivnlivcHi 

(q) coolont flow equations; 

(b) heat transfer equations; 

(c) neutron dynamics equations; 

(d) equations for external effects. 

I hese equations with these or those suijpositions are solved, as a rule, hy means of 
computers. For the quicUperating part of the process may Im used the method, accounting 
for the distribution of parameters (along the reactor hciglit). The hydrolic port of the problem 
(the change of coolant circulation forono-phase liqaid) is usually solved independenUy from 
the thermal one. The neutron dynamics ot the intermittent insertion of large negative 

reactivities may be also separated from the thermal ones for the quick-operating part of the 
process. 

For the pur|)oses of preliminary designing, as well as for a Imttcr understanding of 
thermophysical and hydrodynamic processes proceeding in reactors, the analytical dependences 
are also useful which give the possibility to define in a general form the characteristic 
criterion of the process and approximately find the interconnection required for limiUng the 
transient process within the tolerable limits of parameters changes. 

Below, the main res lits of the study are outlined and some recommendations are given which 
should be followed, to our mind, during designing the reactors. 

The more detailed data are supposed to be published later. 

2. SYSTEM OF EQUATIONS FOR EXPRESSING CHANGE OF COOLANT CIRCULATION 

IN REACTOR CIRCUIT 

The coolant flow equation can be easily obtained from the mechanics equation: 


(mw) 


and it has the following form: 




Using the appropriate approximations for the foroes of friction, pressure, gravity and local 
resistances acting in element dx of the circuit , from formula (2.2) it is possible 

to obtain the equation for time change of non-condensable, coolant circulation in the closed 
circuit in the following form: 

“ A Pmov - A Pres ? (2.3) 


where: 
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(the summation is carried out by oircniil sectors with mean parameters). 

Adding the well-known cciiiation for chun^a? of an^^ular velocity of the puin|) rotor rotation to 
formula (2.3) we shall obtain the followiafi; system of equations: 


at 




i - M M 

jjp “ ‘'hnov '^res * 


It is convenient to write the system of equations (2.5) in a dimensionless form: 


„ , d /G , ^^mov 


rn>och^(a 


d ((o\ l ^mov . ^res 

ir~ ‘IT 


where-the values: 


*^^0 /I V \ ^21^0 . 


^mech =■ 


I^uq 2nif 


have the sense of characteristic time constants of the transient process. 

Value r|^yj(the hydrolic time constant of the coolant circulation change) characterizes a ratio 
.of the kinetic energy accumulated by the coolant to the power required for overcoming the 
coolant flow resistance* 

Value rmech mechanical time constant of a change of the rotor rotation angular velocity) 
characterizes a ratio of the kinetic energy accumulated in the pump rotating parts to the power 
in the pump engine shaft. 

From equationg (2.2) a similar system of equations for the reactor with the branched 
primary circuit may be easily obtained* In particular, for a reactor with n loops in the i-th 
loop of which Pi pumps (provided for the simplicity with common coupling points) can operate 
in parallel, we shall obtain: 


f., ^*’^ik 


Mmovjk ^resj,^ 


^ Pi n 

^1*2 Gji^ (k«l,2,...Pi) ! Gp-S Gj (iipl,2f..n) . 

k«i 1-1 


Physically, the apparent conditions of the operating consistency of pumps in transient 
processes are as follows: 

382 
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(n) prossure oqunlity in tliu tniiiHient |>m(;oHs in conmion points of ilio fusion of different 
coolant flows; 

(b) continuity of all tlie flow nitos nml nn^i^ular vc^locitios. 

l*or illustrating [)ur|)oscs lot us consider tho comlition of the openitin^ consistency of pumps 
in the i-tli loop, which is easily obtained from fommln (2.9) and has tho followin^r form: 

movik -M’resjij -t//ik -‘ tl t -icl(«m(t). (2.10) 

dt 

After disconnecting? of the ik-th pump, condition (2.10) is automatically maintained at the 
expense of the ''kinetic force" flow rate till value has reached 0. When starting 

the pump, condition (2.10) may be satisfied only in case the started pump develops a certain 
pressure -head, Up to this moment the* pump may be considered to be operated in the non-delivery 
conditions. According to equation (2.10) the pump pressure head should not be less than a 
difference of the pressures in the inlet sector. 

The conditions similar to condition (2.10) will take place for the joint operation of pumps in 
different loops as well. 

It should be noted that if there is no flap in the pump, the reverse flow of the coolant 
through the stopped pump is possible* In this event the equations should be changed. 

At last, it should be noted that after the pump is deenergized the interia of a jet can become 
an external force relative to the pump wheel. In this event, the pump will be under the 
hydroturbine operating conditions. It is obvious that it will take place in case the following 
equality is obtained: 

^•^movik -0 (2.J1) 

(in a general case* at (Uj|j,4 0). 

From this moment on the pump wheel is transferred from a motive force into a resistance; 
this concept should be put or taken into account in a system of equations, type (2.9). 

Taking into account the considerations outl.ned above, the solution of the system of 
equations (2.9) may be provided by a computer. 

The analytic solutions may be obtained by using a series of simplified suppositions. Below, 
these solutions are given for the case' of interest for the circulating pumps. 

3. CIRCULATING PUMP DEENERGIZING OR BREAKING 

Under the term ’’pump deenergizing” we shall understand the deenergizing of feeders in 
case'the pump is provided with an electric drive. 

Using the well-known statistical dependences it is possible to approiimrtely put down: 



where value expressing the. pump transconductance can be obtained from the pump 

calculations and the experimental data. 
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The circuit cliuracteristics usually approach a square value of the flow rate, i.e*: 


Ap .((i )2 

(‘X2) 

Moment of resistance in the pump is 


''^njK-Mliyd+Mfr ) 

(d.,3) 

where: 

Ml, yd =, ; 

(.'{.4) 

Mfr . 

(3.5) 

For analysis it is convenient to use also the following correlation; 


Mfr -Mfrx a . fi ^ 
a»o Uq 

(3.6) 

which gives the dependence within the range of ^ ^ quite near to expression (3.5) 

Go 

(value Mfj. within the range of £ »S!L is unimportant). 

Gq 

It is considered in the first approximation that from the moment of the pump deenergizing 
we have: 

^mov~0' 

(3.7) 

The effect of M^Qy on the process will be evaluated separately. 

The system of equations (2.6) together with the above-mentioned dependences will be 
written as follows: 

'•hyd 

(3.8) 

„0 l^O 


wherei r- 

x-H. ; 

Go 

(3.9) 

6>0 

(3.1-0) 


From equations (3«8) we 'shall approximately obtain: 

1 


X = 


wherei 
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^circ. (the. time constant of coolant circulation drop) seems to be in the most simple 
case only a sum of the preliminary dcteriuineil partiol time constants. 

When 


^c“Np“Npjj5 


(3. 13) 


(that aiay take place in an ideal case, i.o. ntr)jj*»jvoI "D, 
then: 

r.u. 2(Kc+Rp) 

i.e. value rdrc expresses n ratio of the. kinetic energy nccuiimlnted in the circuit (by the 
coolant and the pump) to the power for coolant pumping. 

The. greater the. accuracy of expression (3.11), the belter ihe corrclulion is satisfied: 


i!lZi ..med, 


(3.15) 


Condition (3.15) is not rigorous. It may be shown that the solution for equations (3.8) 

within the range. of rhyd-v-miech and ,^o-0, corresponding to the worst accuracy of approximation 
for equation (3.11), has the. following form: 

x^y]/l -2lny 

j\—2\wy‘ 


r “l/F' 


or, the solution for the nearest case is: 


2 ( 1 + 27 “) . 

5 r-rmech at — ^ - 1 ; 

(1+21)2+1 ’’aiech ^ 


“ ^mech 


( 1 + 21)2 + 1 


Vd 


. Hi yd 
at — ^ = 


The comparison of dependences (3.11), (3.16), (3.17), made in dimentiotiless coordinates 
^ ^ shows that withifi the range of T^ydrvrfnech expression (3 J 1) gives a considerable 
relative error, especially within the range of small flow rates. However, .the absolute error which 
the designeris Lnterested in does not exceed~0.05 in x parts and this is quite suitable for 
practical calculations. 

Using dependences (3.11) out of the system of equations (3.8) it is possible to obtain an 
expression for y, .which allows to determine the moment of the pump transferring into the 
hydroturbine operating conditions. Below, the following limiting cases are of great interest: 

(a) '■fflech>^'’hyd . 
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In this cnsc (practically, independent of value (4^,) we shall obtain; 



— L 


CiJIl) 


mecti 


i.Gi, the process is determined by the intoria of rotatinf^ parts of the pump. Tliis is the ranj^o of 
proportionality between thc pump rotor anf^ular velocity variation and the coolant flow rate. In 
this case itia possible to employ relation x«y used by many authors (for example^ See [2]). 

When employed in the system of type (2.9) this relation can considerably reduce calculations and 
in some cases will help to obtain approximate analytical dependences. 

^mech»^hyd desirable to determine the electromaf^nctic slowlne:-dowti of a rotor of 
the centrifugal pump with a closed-type asynchronous motor employed in high-pressum water- 
moderated water cooled reactors. After the high power pump is deenergized a slo wing-down of 
its rotor occurs due to the interaction between the lowering rotating magnetic flux coupled with ilio 
rotor and the eddy currents induced by this flux in the pump stator material. Taking into account 
that the eddy current losses power (Np) is proportional to the square of the number of full 
cycles of magnetic polarity reversal (in our case — to the square of the pumu rotor angular 

_ngt 

velocity) and to the square of the greatest value of magnetic induction (e )2, it is possible 
to obtain the following expression: 

^2t 


/el.mag 




circ 


X = 


/ ’’f 

el j t 


^el.mag 
circ 
mag 


*2l 


circ ^el.mag 


where: 


Np 

■IV 


'’el.mag "■ 




(3.19) 

(3.20) 

(3.21) 


(Lf^ and Rr — the inductance and the resistance of the rotor winding, correspondingly). 
For the most practical cases value «1. 


circ 


Having considered that in this case the main part of the intergrand is valueJL, we shall obtain 


from formula (3.19): 


-2t. 

1+, ,^l»tnag /el.mag 


X sss 


(3.22) 
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Thu=, . dr.,, „t .h„ rl,m p|„„ to ,l.o 

additionnl Hlowinj?-<i()wii of the rotor cIooh not oKcei'd, intliis case, the value of ,.r .1 m 
that is iiHimlly not inorc; than 0.07 in piirtH from x. *' 

hxproM.sion (.1.22) indicutc.s that them! effoiit.s decitni.sc with un increase of value r ,• rs- r , , 

At value r ,, = (2_;|) sec. they may iKMieglectecl. (Isunlly, valuer, „ech^ I h«^c. Addil^mal'""^ '' 
pump nywheel ma.s.ses may con.sideralily increase value 

", 'moch»1„,l ..r th., „l„„, i, ir|Lr»illc Jete,mi.e 

. rdottvedm,, „r ll,c Ito, ™,c i„ „oc„.l.„co wltl, ,l„. 
rclative rotor un^riilar velocity. 

'■niech»»tiyd of the. test bench it is |)o.ssibIe to make mcasurctnents on an unnatural bench 
t . should he homo in mind that according to formula (2.8) for hot and cool systems there is the ’ 

followxn«a„,mxin,ate.correlation: rJe,,. ..r^,, (:i.23), i.e., a dmp of the coolant 

flow rate in the hot circuit will iKi more ^jradual. 

It is noteworthy that for thcreactors with r,„ech»rhyd » niechanical breakinBor jamming of 
t e pump rotor may be dangerous. In thus case an accidental drop of the. coolant flowrate 
may occur which is detemined by low value r^yj. b'or these reactors it is desirable to have an 
emergency pump constantly connected (in parallel) to the cirouit as a precaution measure* 
'■hyd»rmech 

In this case we shall obtain: 


~ 1 
x=: 

1+-L 


) t<to 


(3.24) 


(^ IS the moment of pump transferring into the hydroturbine operatinR conditions). It may be 
shown that in this case value y tends rapidly to the following expression: 


t'O 

) 0<<^o<l 7 

'’hyd 


(3.25) 


i.«., .cco*B 1 „ ( 2 . 1 1 , a„d ( 3 .,, ,he p™p h,dn>,«,bi„e.„pep„i„B 

OTd.uons. ft m, be assumed with a mapja that rrom thetnoBent at (.amp deeaargising the 
following correlation takes place: 


where: 


... APwheel 


Value is of the order of 0.3 - 0.6 and may be obtained daring the pump testing. 
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(It should ho iiotnd lliat llio omployimMit of ndiitioii x-y for simplifyinf; the culrulalions is 

physically not worth whilo in this coao. Itut it ounnoi ItMid to an ermr if — «l' 

AP„ 

It is explained liy the fact that at niy(|»r||ni(;h the contriluition from tlie puni[) I'otating |)nrts 
is ncKliKiblc and the upproxiiiiation nalure of value y for detenninint' x is inessential). 

In this case the cx|)t?rimantnl determination of a drop of the coolant flow rate seems to Ije 
possible only in a real circuit. It should be noted that if tests are carried out in a cool circuit 
and the ratiiiR is self model then it may be considered that according? to the eqiinlioa (2.7): 


fl' .wrtJO 

hyd hyd ’ 


(3.2«) 


i.Bi, a drop die. coolant flow rate will be almost the same both for the hot and cool systems. 

The analysis of home big ssationary reactor plants (for example, of the WWPH type) shows 
that the. determining time constant for the considered process is value r^jyj. 

Qualitatively, .it is possible to write: 

7 ■ 2 

w (A . if ^ (:i.29) 

"'ar , "1 — >' 

“hyd 

The average. speed of coolaat in the circuit (wq) and the. quantity of local resistances in the 
reactor per a length unib(if(b)' slightly change for various water-moderated water-cooled 
reactors* 

In this case, usually 


Therefore, for the reactors with large liydrolic diameter (d|,yj) value rhyd is relatively great. 

For the small-size reactors (for example, of tlu-*.icc-hrcaker "LclNIN” reactor type) the. 
process is mainly determined by value r^ech* analysis of value rmech is difficult since 
dependence I(Np, a>o) is strictly constructive. When designing the. pump it is expedient to put 
forward the requirement to have, if possible, value rmech great enougli being ?(1 A 1.5) sec., 
which is possible. to provide. 

At last it should be noted that values r[,yj and rmech order of (0.2 t? 1) sec. 

(The latter is obtained when the pump has no attached flywheel masses, for example due to 
the external pump engine for setting a low pressure). 

4. TEMPEHATURE PERFORMANCE OF REACTOR IN TRANSIENT CONDITIONS 

Equations of the reactor thermal dynamics for slow processes or the processes comparable 
with the time of establishing the thermal equilibrium in a reactor may be obtained by reducing 
the -corresponding equations in partial derivatives to the mean values. 

In the raost simple for-^. they are reduced to regular balance equations, for example* of the 
following form: 

Cu -^-(qoF)^ -(kF)u(Tu-Tcan); (4 1) 

dt “O 


-^9 ^ 
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dT. 


Lean ^ can . Wkr )can(^cQn~f ’c)/ 


p die 


^ ^ ^ out-'nn)> 


wherej 

^out)* (^*2) 

Tin -fU) (4.3) 

and'Value Tc may be given proceeding from the physical considerations. 

We sliall not consider tlie-equations of delay and heat exchange in communication lines 
and steam generators because for the most interesting quick p»art of' the process in the given 
problem the -operating conditions of reactor and steam generators are [)ractically separated in 
timei 

The system of type (4,1) together with the above-mentioned equations of type (2.9) for 
changing the coolant circulation and the known neutron dynamics equations give tlie 
possibility to express (in the point approximation) the transient behaviour of the reactor. 

To find just the thermal time constants of the transient processes in the reactor, let us first 
consider transient thermal process when the emergency protection system is tripped at 
a constant flow rate of the coolant. This condition has much in common with the transient flow 
rate problem and is of particular interest. At n jump introduction of the great negative 
reactivity (pem.prot the emergency protection system, the temperature reactivity 

effects for the quick part of the process (i.e. at T|^ Const) may be neglected. 

Then, at one group of delayed neutrons we shall approximately obtain: 

t 


c 


J3. 


^del 


Vq ^+Pem,prot 


(4.4) 


It can be shown that the effect of fuel elements canning on the character of transient processes in 
the watet^moderated watercooled reactor is small (Ccan.«Cui Ccan«Cc). The value of canning 
heat capacity can be effectively considered in other constants (for example in Cy). 

Then the system of equations (4.1) including formula (4.4) can be reduced to theTollowing 
form: 


'^pas 


'^u 


£o 

dt2 


where 1 




Tn-Ti 


in 


rp O 

^ c*" Mn 




at Const«l ^ 


• Const— 




^del 


em.prot 


^(4.5) 


0*^6 r 


Tu-Tc 
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^u Qu 


— T 


at Const-1— 


as 

^del 




N. 


(4.6) 


(4.7) 


(4.8) 
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Lc + 


“ .s|i — " V" * 

V I V tQ 

Value r^, - the time constant ofcliiuigin^ tlu* fuel elcMiient (unit) l(un|)<?mtun‘ - rlmrnct(‘ir/.(*H 
Q ratio of the excessive heat {i\s compared to the coolant heat) accuinulatetl in the unit to the 
reactor power. 

Value rtf — the transport or cunvection time constant of the reactor — charactcrize« a ratio 
of the heat generated by the fuel elements and coo kmt when they are cooled from the coolant 
average temperature to its temperature at the reactor inlet to the reactor power. 

In the presence of the equation: 

'’tr rpa s * ( I + (4. 1 1 ) 

i. g 

we may say that value also characterizes the delay in the convectionai heat transfer (along 
the reactor lengthjby the coolant due to the terminal heat capacity of the unit, i.e., it is the 
time for a heat pulse to expand in the reactor. It should be also noted that in the energy 
content of the reactor time constants considered herc^ value 

^tr Qnc 

(4.12) 

is the time constant of establishing the thermal equilibrium in the reactor and characterizes 
a ratio of the heat stored in the reactor to the reactor power. 

Equation (4.5) may be easily solved. 

We shall be interested in some limiting cases which are of a practical interest: 

(a) Then from equation (4.5) we shall approximately obtain 


^“^Pem.c 


Pem.prot 

^■^^em.prot 


i.e*, the transient process is determined by the time of heat wave expansion in the reactor. 
Physically, it means that the excessive heat in the fuel elements as compared to the coolant 
is neglected. In this case the fuel element temperature is close to the coolant temperature and 
it is possible to speak about cooling (heating) the reactor as a whole. 

(b) ru»rtp. 

In this case from equation (4.5) we shall obtain: 

0+Pem.prot ^'*"^em.prot (4.14 

i.q*/ the transient process is determined by the thermal inertia of fuel elements, and a 
relativc'drop of the-coolant temperature approximately coincides with a relative drop of the 
heat flow,. Physically, it means that in this case the coolant heat capacity may be neglected. 
382 
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(i follows fmin tht> oxnniploH oonsith'nid abovi* tluit tli(' pump stoppn^o coiulitions are mart? 
(InngorouH for the reactor with iho fuel (liaMioiits of a lii^h ihonnal polentiul, 

(this, for oxfuaplo^ occurs with the fuel (‘loinoats inudo of .sintcrtHl uranium dioxide). 

It is possible to write* that: 


r 


^canl’u ^cl'a 


whercj 


(A — the coefficient of unit shape — may be obtained by solviiifi; the equation of heat 
conductivity in steady conditions). 

For the -medium size units using uranium dioxides (for example, for the ice-breaker ^*LKNIN** 
reactors), the first two expressions for will be the main ones (due to the small heat- 
conductivity of uranium dioxidc-and the gns gap). 2 

For the large size uranium dioxide units value is . Usually it is t ^^{2 f- 5) see. 

for the uranium dioxide fuel elements. 

It is easy to show that when satisfying the condition: 

If- 4 (4.17) 

''circ ^■*'^ern.prot 

the coolant temperature in the reactor will not rise just after the pump deenergizing with 
a simultaneous triggering of the emergency protection system. 

Here: 2 

..r f^U 


Physically, condition (4.17) is quite simple* it requires that the initial drop of the heat flow 
should be more sharp than that of the coolant flow rate* For a deep drop of the power 
(^em.prot»0) condition (4.17) will be the following: 

^circ* (4.19) 

As mentioned above* value rc[j.c'--(l-1.5) sec.for the reactors where the pump is not provided 
with additional flywheels. Therefore condition (4.19) for the uranium dioxide units may be not 
completely satisfied, i.e., in this case a considerable-rise of the coolant temperature is to 
be expected especially when there is delay between the signal arrival to the emergency 
protection system and the-beginning of its efficient effect on the reactor power change 

(fem.prot equals some tenths of a second). 

An approximate calculation of the parameters distribution (with respect to the reactor height) 
can be-made in the-following way. 

Having accepted that in the -heat exchange-equation for the moving liquid 


-12 - 




tlic Vtiriiiblfs in the expreHsion for ([j, urc He|iiirat(ui, it in (K).s.sible to write its Moliilion as 
follciwa. 

(u) l•’o^ the coolant |)urticles l)oia(< iit the initial moment of tranaiont proceas (t-O) inside the 
rone tins core (t,’j^0): 


;(x,t)-iHt(x~fw(t")dty 0)+/qj(x')J^ . , 


where: 


x-/w(i"),|t’ 


x'-x-fw(t")dt" 


(4.21) 


(4.22) 


x-fw(t")(lt" 

ist(>^“fw(t*‘)dt7.0)-i . f.|-.‘(x'‘)dx»‘ . (4.23) 

b 0 I 

(b) Por the coolant particles beinp: at the initial moment of transient process outside the 
reacting core (t^^O): 

ic(Xjt)«iin(o,ty + fq (x^) . 

C" in' » O' I I] ' ' No G(t') ( 4 . 24 ) 

(Moment t^ of a coolant particle entry in the reacting: core at the point of interest (: , i) is 
determined from expression (4.22) at x^=0). 

The relative value of the heat flow transferred to the coolant 


N ^ (kh )can.(T"can""Tc) 

^*^(kF)o (TO _Tg) 
can can ^ 


(t.25) 


can be obtained from the problem in point approximation. 

It should be noted that for the fuel elements with a low heat potential (r^jc:<r[|.) the essential 
simplification of the problem is possible, namely: in dependences (4.21) — (4,24) it is 
adequate to supersede the coolant speed with the speed of heat wave expansion: 

’ (‘i-26) 

In this case a relative drop of the neutron flow may be approximately used instead of the 
unknown heat flow. 

In accordance with the plan considered above, this problem may be solved both in the point 
approximation and by calculating the parameters distribution. But these solutions are bulky. 
For preliminary estimatipns (marginly) the point model may be -employed in its most rough form, 
assuming that: 


^out “ Hn + 


N 


in + s* 


(4.27) 


It can be shown that for the elements with ru»rtr for which tlie given problem is actual, 
a drop of the heat flow at a decrease of the coolant flow rate up to the flow rate of a small 
(steady) circulation practically coincides with that at the constant flow rate. It is connected 
with the fact that a drop of the heat flow due to the rise of the thermal resistance between the 
coolant and the unit is small in the given case. 
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111 view Ilf tlie nliovc-iiientioiied fiiiitoi’H tliu coiiiiection of tlio proccsH clinriirtoriHlic 
piirniTUJtorH cun bo easily olitained al wliich tlio relative lemiioratiire (or the eiitlmlpy) of the 
noolmit will not exceed tlio Kiven value flu,,,, for example, in the followintr form: 


where: 


— r'circ 


Mim “ 'ill 


•'’om.prot 


Kl^ peiii.prot 


(4.2li) 


im - )l» 

/ I'ui. prot/ff 


Value fl|j,„ is apparently to be taken from the no-boiliiiR condition of the coolant in the 
reactor. 

Value Pem.prot/fl taken from the temperature effect overlapping condition at the reactor 
cooling up to valuc'^Tipis from 2 to 4. In this range it does not limit the process considered. 
The effect of Pem.prot/l^ decreases with the increase of value tq because in this case the 
process is determined to a greater extent by the excessive heat of the fuel elements and not 
by the internal heat release. 

For the fuel elements with r^«Tif practically there is no problem of emergency heat 
removal after stoppage of the pumps. 

In this case condition (4.17) for the absence of the coolant temperature surge is satisfied. 
For the reactors of this type there are mainly two moments: 

(a) Value fem. prof 

With some margin it may be considered that for limiting the coolant temperature' overrising 
by value for time*^rem. prot approximately rated power of the reactor) it is 
necessary that: 




(4.31) 


For the intensive power water-moderated water-cooled reactors value taken from the 
no-boiling conditions of the coolant, is equal to 2-3 and condition (4.31) may be easily 
satisfied (if there is no breaking in the pump). 

It should be also noted that it is desirable to have: 

^erntprot*^ '■pa^ ’ (4.32) 

In this case the integral of the heat accumulated by a coolant particle -for time^rejn. prot 

will be small. 

(b) To provide the given value in the process end it is necessary (marginly) that: 




(4.33) 


illll 

prot 

At last the approximation of type (4.27) is noted to be essentially improved by selecting the 
physically approved effective centre of heat removal in the reactor and to be reduced, for 
example, to the -following form: 
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The. express ion of type (4.34) iiuiy give tlio close coincidoiice with the problem of distributed 


parometera (along the reactor height). 

5. IlhlAT EXCHANGE CRISIS IN REACTOR IN TRANSICNI’ PROCESS 

Using the approximate static dependences for determining^ the critical heat loads it is pos- 
sible to evaliinte the margin on a heat exchange crisis in the transient conditions, f-or 
example^ using (at a constant pressure) the depondeuce presented in f3]' 

1 1 

qcr = Const (wy) ^ x (Tg - Tc)^ ^ 

it can be shown that for the considered problem the margin on the critical load in the 
transient process will immediately decrease, if; 

!iL 5" 2 prot 
'"circ prot 

in this case the value of small (steady) circulation can determine only the minimum to which 
this margin will drop. 

It follows from this that when satisfying correlation (5.2), it should be given that: 

”q„. (5.3) 

'^trans 

i.e. it is necessary to have a dynamic margin already in the stable operating conditions. 

It should be noted that the transition into boiling conditions will require the further increase 
of this margin. 

From the previously given data it is also clear that correlation (5.2) will, for example, 
take place for the fuel elements made. of uranium dioxide. The transient conditions may be 
preliminarily estimated (marginly) in the point approximation. 

Let us assume that: 

(a) Position of the dangerous point in the reactor (x^an) is not changed in the. transient 
conditions. 

(b) Minimum margin on the critical load (nqcr) occurs at the moment of transition entry 
into the stable (small) circulation. 

(These assumptions correlate quite well with more accurate calculations including the 
parameters distribution along the reactor). In this case,formula (5.1) including, the following 
parameters correlation may be obtained which does not result in (marginly) the. heat 
exchange . crisis in the transient conditions connected with stoppage of the pumps and 
operation of the emergency protection system: 
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^s(X(Jna)- 


wherci 


^iTm ^ Trans 


^s(xau„)- ^ 


r«^T 


o 

In 


( 5 . 4 ) 


(5.5) 


JV“(*dan)-T|{, 

value gein expression (5.4) being taken from equntioii (4.28). To cleteriiiiiie q^r the dependences 
of type (5.4) can be obtained with approximations different from equation (5.1), If 
estimations should be made frequently, dependences (4.28) and (5.4) nmy be shown graphically 
in the form of nomof^rams. 

If the condition of type (5.4) in not satisfied, the heat exchange crisis may occur in spite 
of the chain reaction cessation. In this case the high heat load in the reactor is very 
specific because.it results from the high temperature potential of the fuel elements. At the first 
moments after a steam blanket is formed on the unit surface* the temperatures are equalized 
along the core section and the element canning. The further cooling (heating) of the unit as a 
whole depends on the correlation between the fuel element temperature drive, its internal heat 
release and the heat removal after the heat exchange crisis. 

An approximate. analytical solution of the problem may he obtained from the equation 

system of type (4.1) on the. assumption tlial tlie coolant temperature is constant while at the 

moment of heat exchange crisis the heat-removal coefficient with regard to the coolant 
aiscontinuously changes up to some constant value • 

Then for the most dangerous case (Ccan«Cu) froiS the point of view of the unit melting, we 
shall obtain for not too small values of t: 


'can 


r}T 


where! 


^can = 


"^81 

T'can ““ 


•)e 


,bl 


int 




can = 


n-T° 


u 

, 

'1 'IK 


«r 

'2 


r«r 

u 


TO 

can 

J*0 

can c 

♦ 



; r2 = 

a 

^®can*Fu 


^"^can* b'u 


C-u 

■3P7 


’25 


Cu 


‘can^u 


ei 


jffr 


=ri + r2 


362 


n*=. 




_ t’ 


(5.6) 

(5.7) 

(5.8) 

(5.9); (5.’n) 

(5.11) 

(5.12) 


4m 0 + Pem.prot 
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Vcr ~- 


'•] + '2 


^2 

■+ 


81 


n-^r2 

(* is iriK' for tlie lirno moment just preccuKlin^ llit^ Iwnil (^xc:li{ui^() crisis), 


(5.1:0 

(5.14) 


l"rom equation (5.6) we shall obtniu that the initial sureje of value 6^cnu f^ives the following 
value: • 


^can (^) ^ 


’’2 

r2(rl+ M) 


(5.15) 


At the complete cessation of heat remova I -> <«•) from equation (5.15) we shall obtain: 

'Wf-Tc 

■'can 


2 


^^can M —0^ 


r2 


* rnn"” * C 


(5.16) 


which means that ^ Tti, i*e*, at the period of temperature equalization along the unit 
the canning temperature will not exceed the average core temperature occurred at the moment of 
heat exchange crisis; physically, this is quite obvious, 

A further behaviour of is determined, as it is seen from equation (5.6) by the relation: 

(5.17) 




For example, to lower the canning temperature after a short surge for the most dangerous 
(theoretically) case — the heat exchange crisis at the moment of scram system triggering it is 
necessary that: 






prot 
bl 


(S.lfl) 


'rhe analysis of real values of shows that for the elements with ru»rfr (fcr example, 
elements made of uranium dioxide) condition (5.18) may be satisfied. It means that canning 
temperature in the transient conditions does not exceed (marginly) the average core 
temperature occurred in the rating. Therefore, if 


•pmelt 


(5.19) 


9 

can u 

then, the melting of canning for the fuel elements of this type may not occur. Condition (5.19), 
as a rule* is satisfied on the uranium dioxide fuel elemets for the high-melting materials of the 
canning (steel, zirconium alloy) at high enough heat loads. For the fuel elements of this type 
with the aluminium alloy lower heat loads are essential. 

If ri«r2 then, after the temperatures equalization along the unit, a further increase of 
value 


( 5 . 20 ) 


is possible* In this case it is necessary that: 

/elt nfnt 

f/can * 

%1 

The similar conditions should have place for the core temperature of the fuel element as well. 
Thus, in case of the heat exchange crisis on subcooling water the conditions may occur which 
do not lead to melting of the-fuel element canning or core* 
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The detailed nnolysis when designing the reactor may hel[) to deteriiiliic the permissible linilu? 
of short-time operation of the reactor within the ranRe of increased temperatureH of tlie fuel 
elements. 


MAIN SYMBOLS INDICES: 

G - weip;ht flow rate of coolant mov - moving 

0 ) - anjj^ular velocity of rotation mech — mechanical 

I - moment of inertia vol ~ volume 

M - moment of rotation fi* ““ friction 

m — mass c ~ coolant 

AP - pressure drop , P ““ 

^ ^ el. mag - electromagnetic 

- power I ^ refers to length term 

E - k,.et.c s. - stead, 

C - heat ammount 

,,, , , s - saturation 

W - coolant speed 

fi . r 1 u - unit 

S - passage cross-section for coolant „ , 

F - heat exchange surface *“ houcault 

d- diameter can - canning 

n - perimeter , 

UL-leagth cr‘-“c°risir 

A - local resistance coefficient . , . . , 

. bl ~ blanket 

rj - efficiency . , , 

, ^ int - internal 

n - neutron flux 

p - reactivity ^ 

i8 ~ delayed-neutron fraction me t - me ting 

i- enthalpy pas - passage 

T - temperature “ limiting 

C - heat capacity ““ ^^^^y 

y - specific gravity ““ circulation 

a - heat-removal coefficient element 

k - heat-transfer coefficient ~ 

a - thermal diffusivity 

q- heat flow P^"^ - pumping 

g - gravity acceleration protection 

e - 2.71828 ...... constant sp - specific 

r- time constant ac - accumulated 

. . . eff - effective 

t - time 

trans - transient 
dan - dangerous 
rat - initial conditions (rating) 
r - reactor 
res - resistance* 

REFERENCES hyd - hydraulic 

7. L^F .Epstein^ Journal of Nuclear Energy, Part B, V. 3, 1962 
2, D.LBlokkintzev, M.EJiinashin, U.A.Scrgeev^ '^Atomnaya Energiya*% No.l, 1956. 

5. B, A, Zenkevich^ V,LSubbotin^ '’Atomnaya Energiya*% No.S, 1957. 

4. CuM^Boyd, J .r.yRJLPosseij B^B^Cardtvcll, /.r.» Nuclear Science and Engineering, V, 9, No.4, 1961. 
So K.Pfleiderer^ Mashgiz, W*, 1961. 

6. B^V^KitaeVf l.F^Grevtzev^ Sudpromizdat, 1960. 
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time constant 


t - time 
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